Effects Associated with the Thermal 
Response of the Tl Telephone Transmitter 

By C. A. FRITSCH 

(Manuscript received April 16, 1968) 

The thermal response of the Tl transmitter, when excited by a bias 
current, is analytically obtained. The thermal expansions which produce 
a decrease in electrical resistance are described. Good agreement with 
experimentally measured temperature rises and displacements is demon- 
strated for the early time transients during which the resistance change 
occurs. The displacements which produce "thermal packing" are found 
to be a stronger function of the thermal expansion coefficient of the dome 
electrode than of any other part of the transmitter assembly. 

I. INTRODUCTION 

The design of an effective telephone transmitter requires that an 
acoustic signal (voice) be efficiently converted to an electrical output. 
One method of performing this function is to use a moving electrode 
(attached to a diaphragm) as one wall of a chamber containing 
granular carbon. Thus, if a dc bias current is impressed on the trans- 
mitter the resulting variations can be used to operate remote telephone 
apparatus. 

An important factor in the design of such a transmitter is the con- 
trol of thermal effects, not only caused by variations in ambient 
conditions but also arising from the heat generated each time the 
telephone set is connected to the line. These thermal effects result 
in dimensional changes in the transmitter body, resistance change 
resulting from temperature rise of the granular carbon itself, and the 
resistance change caused by thermal expansion of the carbon granules 
when they are heated by the biasing current. The total effect, which 
produces a loss of carbon transmitter efficiency, is referred to as 
"thermal packing." 

To assay the relative importance of the various thermal effects it 
is first useful to ascertain whether the carbon granules should be 
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treated individually or as a continuum. It has been observed that the 
contact points resulting from the granule asperities cause a "bunch- 
ing up" of the current. 1 Consequently, for ac heating in the kilohertz 
range, the high local rates of heat generation produce larger tempera- 
ture changes at the points of contact with a proportionately large 
drop in electrical resistance. 1 However, it can be readily shown (see 
Appendix A) that the carbon particles are small enough and their 
thermal diffusivity is high enough so that any local temperature 
changes are virtually diffused in about one millisecond. Since the 
resistance changes in the Bell System Tl transmitter* have been 
observed to take place in about one second, 2 then the local heating 
caused by the asperities and the accompanying resistance change can 
be neglected and the granular carbon can be treated as a continuum. 
To determine the effects of thermal expansion in both the micro- 
phone body and the granular carbon itself, the temperature distribu- 
tion throughout the transmitter must be known. Thus, the work 
reported here consists of a first-approximation type of analysis to 
determine the temperature distribution in the carbon chamber of the 
Tl transmitter. The associated thermal effects are then considered 
with the hope of providing a better insight into what has been ex- 
perimentally observed. 

II. ANALYSIS 

2.1 Thermal Analysis of the Carbon Chamber 

The telephone transmitter design of interest here is shown in Fig. 1. 
We wish to describe the transient thermal response of this transmitter 
as the above mentioned dc current is turned on. Geometrically, the 
carbon chamber consists of a dome (moving) electrode connected 
with a conical back electrode by a flexible nonconducting chamber 
closure. If we assume that the surfaces of constant voltage and con- 
stant temperature within the carbon aggregate are hemispherical then 
the walls of the carbon chamber can be considered as two concentric 
hemispheres (see Fig. 2) . 

If we also assume for the moment that the relatively heavy back 
electrode is held at the initial and ambient temperature T , then, the 
thermal response of carbon and its chamber can be conveniently 



* Made by Western Electric Co., the manufacturing and supply unit of the 
Bell System, and available only to the Bell System. 
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Fig. 1 — Cross section of Tl transmitter, 
written in terms of a temperature excess 6 defined by 

6 = T - T . (1) 

If we further assume that the temperature excursions are small 
enough so that no significant variations in the physical properties 
take place then the energy-balance equation is given by 



pc a7 = kv * 6 + q ' 



(2) 



where 



k is the effective thermal conductivity of the carbon 

pc is the heat capacity of the carbon 

q'" is the rate of heat generation per unit volume. 

Because of the angular symmetry of the boundary conditions and 
the heat generation, gradients in the directions of the angular co- 
ordinates can be neglected. Thus, on dividing through by the thermal 
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Fig. 2 — Coordinate system for carbon chamber analysis, 
conductivity k, equation (2) in spherical coordinates becomes 



136 

k dt 
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dd 
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k 



(3) 



where k is the thermal diffusivity. 

In general, the heat generation per unit volume is given by the 
product of the square of the current flow per unit area times the 
electrical resistivity. Since here we have only radial flow of current 
7 then 

q'"{r) - (I/2^r 2 ) a p. = i 2 aPe . (4) 

It can be readily shown that the resistivity p e is related to a total 
measured resistance R for a cavity formed by two hemispheres of 
inner radius, a, and outer radius, 6, by 

Pc = 2T[ab/(b - a)]R. (5) 

To specify the boundary conditions on our problem we recall that 
the initial and back electrode temperature are taken to be ambient. 
To approximate the heat lost by the dome electrode we specify a 
certain thermal resistance between the dome and some sink at am- 
bient temperature and that resistance is represented by a coefficient 
h in the so-called "radiation" boundaiy condition. Consequently, on 
using (4) and (5) in (3) the following boundary-value problem can 
be stated. 



I §1 - lF<}?( 2"^V-l 2 ( ab \n 
k dt ~ r\dr V l dr/r ? \b - aF 

6(r, 0) = 

e(b, /) = o 

k dO/dr(a,[l) = h6(a, t). 



(6) 

(7) 

(8) 

(9) 
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where 

= 7 2 22/4ir/fc. (10) 

The above nonhomogcneous problem in spherical coordinates can 
be simplified by the following substitutions. If we define the dimen- 
sionless variables 

r = Kt/(b - a) 2 , (11) 

V = (b - r)/{b - a), (12) 

and 

, ( ,, r) = x&JL + ^L. (1/r) , ( i3) 

then the problem specified by (6) through (10) becomes a homo- 
geneous transient conduction problem in a slab with an initial tem- 
perature distribution: 



dv dp 
dr drj 



(14) 

»(,, 0) = F/(l - F)( x _ ( /_ Fh ) (15) 

v(0, r) = 27(1 - F) ( 16 ) 

dv/dr,(l, t) + i1 ^/ 1 {N Bi + 1M1, t) = (N Bi + 2)/F, (17) 

where 

F = a/b (18a) 

N Bi = fto/t. (18b) 

;V BJ is called the Biot number which characterizes the ratio of the 
rate at which heat is lost at the dome electrode to the rate at which 
heat is conducted to it through the granular carbon. 

The analytical solution of the transformed problem is relatively 
straightforward if the dimensionless temperature v (tj,t) is first divided 
into two functions, one representing the steady-state temperature 
rise, the other corresponding to the transient response. Thus, we set 

v(v,r) = v(v) +»(h,t). (19) 
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The steady-state portion, cp (77) , satisfies 

d 2 <p/dr, 2 = (20) 

-K0) = F/(l - F) (21) 

d<p/dv(l) + (1 ~ F) (N Bi + 1M1) - (N Bi + 2)/F. (22) 

The steady-state solution is then 

The transient portion of the solution results from letting #(»?,t) 
satisfy 

d&/dr = B^/dr? (24) 

*(*, 0) - F/(l - F)[ - _ (1 1 _ F) - ] - *(„) (25) 

#(0, r) = (26) 

30/^(1, r) + (1 ~ S (tf a , + 1)0(1, r) = 0. (27) 

The solution to the problem in #(*7,t) is derived in Appendix B and 
can be stated as follows: 

00 n/\2 _i_ 2\ n\ 

&(v, r) = X) 7 2 ■ " 2 1 ' exp (— X*r) sin X n r? / #(77, 0) sin X„tj dij, 

n-l An T Cl T Ci J 

(28) 
where X„ are the positive roots of 

X„cot X n = -d, (29) 

and the parameter c x is defined as 

c i = {± ^ n (N Bi + l). (30) 

The series solution given above was found to converge too slowly 
for practical evaluation. * However, the problem specified by (14) 
through (17) can be readily solved through the use of finite differ- 
ences. Once v(t),t) was so determined for various values of Nm and 



* The difficulty resulted from the oscillating nature of the integral in (28). At 
small values of t as many as 350 terms were inadequate for convergence. 
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F, then the temperature rise above ambient 8(r,t) could be found 
using (13). 

2.2 Temperature Distribution in the Conical Diaphragm 

In the analysis of the transient response of the carbon chamber, 
the boundary condition associated with the dome electrode has speci- 
fied that the heat loss from every point on the interior surface of 
the dome electrode is proportional to the difference between the elec- 
trode temperature T(a,t) and the ambient temperature T . The con- 
stant of proportionality has been designated as h, which is a measure 
of the impedance of the conical diaphragm to the flow of heat. Thus, 
the heat flux at the junction between the diaphragm and the dome 
electrode is given by* 



q" = h[T(a,t) - T ]-y= -k — 



(3D 



where the coordinate system is shown in Fig. 3. The factor y is defined 




Fig. 3 — Coordinate system for heat flow along the conical diaphragm. 



as the ratio of the surface area of the dome to the cross-sectional area 
of the diaphragm at s = Sj . If the thickness of the diaphragm is 5 
then 7 = a/ 8. 

In actuality, the coefficient h as defined by (31) is not a constant 
since the dT/ds\,„., divided by [T(a, t) — T ] is still a function of 



* It can easily be shown that the heat losses for the dome electrode to the air 
in contact with it are negligible. 
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time, somewhat larger at early times than at steady state. This means 
that the thermal impedance increases from some initial value in propor- 
tion to the loss of the capacitive effects as the transient dies out. Thus, 
the steady-state value of h will be a minimum value. 

A more exact approach to the problem would require the simul- 
taneous determination of the transient response of the diaphragm 
linked to the carbon cavity by a matched heat flux boundary condi- 
tion. However, if the capacitance of the diaphragm is small com- 
pared with the capacitance of the carbon cavity* then the determina- 
tion of h using the steady-state temperature gradient in (31) should 
be at least in the correct order of magnitude. Because of the crude- 
ness of other approximations it was felt that the more exact ap- 
proach was not warranted. 

When we perform a heat balance for steady state axially symmetric 
heat conduction in the conical diaphragm, the following boundary 
value problem can be stated: 

d/ds(s dT/ds) = (32) 

T( Sl ) = T(a, t) (33) 

T(s ) = T . (34) 
The solution is 

T(s) - T(a, t) In («/«,) „ 

T ~T(a,t) In («„/«,)' k ' 

Differentiating (35) and using the results in (31), we then have 
for the steady-state approximation 

8k d cos 12.5° ,„„> 

h = — 2~j — i i s , (oo) 

a In (sM 

where 

s t = a/cos 12.5°. 

2.3 Displacements Resulting from the Thermal Response 

To calculate the displacements of the dome electrode we have 
chosen some rather simplified models for the geometrical configura- 
tions and constraints. 

First, consider the dome electrode itself. If we assume that the 



* The capacitance of the diaphragm is about one tenth the thermal capacitance 
of the carbon chamber for the Tl transmitter. 
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edge in contact with the diaphragm is essentially free it can be 
shown from both physical reasoning and the theory of elasticity 
that for an unconstrained hemisphere with a uniform temperature 
change the change in radius is given by 

Ar = a t d(a, t)-a. (37) 

The calculation of the displacements in the conical diaphragm 
are somewhat more involved because of the temperature gradient 
along the diaphragm. We first assume that the cone angle is small 
enough so that the conical diaphragm can be approximated by a 
disk. Thus, the displacement in the radial direction is given by 3 

u . (x + „) 2l T [T( s ) - T ]s ds + C lS + C 2 /s, (38) 

s J,, 

where 

T(s) is the temperature along the diaphragm as given by (35) 
T is the (ambient) temperature in the stress free condition. 

If we rewrite (35) in the following form: 

m - r = [T(a, t) - iy £|g (39) 

where T(a,t) -T = 6 (a,t), then (38) becomes 
,_ . a,6(a, t) 

w = (1 + " ) tofeAJ 

The two constants can be determined from the boundary conditions 
on the diaphragm. We will assume that the outer edge of the diaphragm 
is fixed so that 

U m 0, S = S . (41) 

If we also consider that the inner edge is free so that the radial stress 
at the location is zero then 3 

du/ds + vu/s - (1 + v)a t d{a, t) =0, 8 = s, . (42) 

Applying these conditions to (40) we find that 

. , a,e(a, t) s . fs\ s 2 t . /s,\ .s s? 

u = (1 + v) ; . . - In I — I - — In (- J - 7 + — 

In {81/80) .2 \*o/ 2s \s / 4 4s 
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-Ms)+fl 



i-ifeJ+l'-C" 



1 - V 



i +e 



(43) 



a + »> 

Since we are primarily interested in the displacement of the dome 
electrode which is fastened to the diaphragm at s = s lt we first eval- 
uate (43) for that location : 



w(si) = 



2s 1 a l 6(a, t) 


hsm-b 


In (s /s,) 


_ -©fe) . 



(44) 



There are two methods by which the displacement along the slant 
height of the conical diaphragm can be exhibited as displacement 
of the dome electrode along the axis of symmetry of the Tl type 
transmitter. As shown in Fig. 4a, if the diaphragm is rigidly fastened 
to the dome then the dome displacement resulting from the thermal 
stresses in the diaphragm are given by 



A^-min = — u sin a, 



(45) 



so that the axial displacement, Ah, is somewhat less than the dis- 
placement along the slant height of the conical diaphragm. Consider- 
ing the other extreme, as illustrated in Fig. 4b, if the fastening is 
loose but the dome electrode is rigid then the edge of the diaphragm 
moves as if up a rigid wall and the axial displacement is amplified 




Fig. 4 — Displacement of dome electrode resulting from thermal stresses in 
the diaphragm. 
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by the sine of the angle. When the growth of the dome electrode, 
a t 6a, is taken into account we have 



A/i mnr = 



u 



a t 6a 

tan a sin a 



(46) 



In the actual case the true axial displacement will be somewhere 
between these two extremes. 



III. ANALYTICAL RESULTS 

The thermal response of the carbon chamber has been determined 
for various values of N B i ranging from to 100 and for two values 
of the shape parameter F. As noted earlier, the conical shaped back 
electrode has been approximated by a hemisphere of radius b. The 
choice of the value of b (and hence F) was somewhat arbitrary. 
Thus, two values were taken for comparison. The character of the 
approximation is illustrated in Fig. 5. It was felt that the smaller 
value, b = 0.44 cm corresponding to F = 0.636, was the better choice. 

The results in terms of the temperatures in the slab, v{-q,r), were 
transformed back to the temperatures in the carbon chamber using 
(13). Consequently, 



ad(r, t) 



F 



1 - (1 - F) v 



[v(v,t) -v( v ,0)]. 



(47) 



zzzzzzzzzzz 






Fig. 5 — Approximation of carbon chamber by hemispherical cavity. Indicated 
radii are in cm. 
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Fig. 6 — Dome electrode maximum temperature rise vs Nm . (^=0.636 for 
standard Tl transmitter.) 

Evaluating (47) at -q = 1, the temperature rise of the dome electrode 
can be expressed as : 



ad{a, t) 
/3 



= [.(1, r) - 1/(1 - F)}. 



(48) 
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Fig. 7 — Temperature distributions at various times for two values of Nm. 
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From the steady-state solution (23) the maximum value of 6{a,t) 
becomes 



(49) 



This function is plotted in Fig. 6 for various values of F. 

The temperature distributions above ambient for F = 0.636 and 
two values of Nm (1 and 10) are given in Fig. 7 for various times 
expressed in the nondimensional time variable t. The value of the 
abscissa denoted as -q = 1 (r/b = 0.636) corresponds to the tempera- 
ture rise of the dome electrode. As expected, the temperature rise of 
the dome electrode is substantially decreased as the Biot number is 
increased. This fact is also demonstrated in Fig. 6. 

The "early time" transient response of the Tl transmitter to the 
"thermal packing" effect will be shown later to be closely linked to 
the transient temperature changes experienced by the dome electrode. 
This transient response is illustrated in Fig. 8 as the ratio of 8(a,t)/ 
6 (a, oo) versus t. The effect of choosing a smaller value of F is also 
shown for various values of N m by the broken lines. It is seen that 
virtually all of the transient temperature change has taken place by 
the time that t = 1. If the thermal time constant r is defined as that 
value of t where the transient has decayed to 1/e of its value, then 
the ratio 8(a,T)/6(a, oo) = 0.632 at t = t . Taking the values of 
t from Fig. 8 then a curve of t vs N m can be drawn for both values 
of F. These are indicated in Fig. 9 by the notation -q = 1. The thermal 
time constant for a position which roughly characterizes the maximum 
temperature within the carbon chamber, that is, rj = 0.7, is also shown. 



IV. EXPERIMENTAL OBSERVATIONS AND APPLICATION 

Experimental observations have indicated that the electrical resis- 
tance of the Tl transmitter is significantly altered when a bias cur- 
rent is imposed. Fischer and Gaudet have attributed the major cause 
of this change to the thermal expansion of transmitter components 
resulting from the joulean heating in the carbon. 2 They tested both 
the standard Tl transmitter and a similar transmitter with a invar 
dome, diaphragm, and back electrode. The low expansion coefficient 
of invar reduced the thermal expansions and hence the invar trans- 
mitter packed less than the Tl transmitter. However, the reduction 
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Fig. 8 — Transient temperature rise of the dome electrode for various Nm 
and two values of F. 

was not as significant as had been expected because the thermal re- 
sistance of the invar parts was greater and consequently higher 
temperatures were experienced on the dome electrode. 
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Fig. 9 — Thermal time constant at two positions as a function of Npi . 
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A detailed description of the thermal effects necessitated measure- 
ments of the dome electrode displacement and temperature. These 
were undertaken by C. Sandahl and P. E. Prettyman. The report has 
not been published but here is Mr. Sandahl's description of some of 
the instrumentation. 

A flexible silicon semiconductor strain gage, 6 mils square by 50 
mils long, was bonded to the diaphragm with epoxy in the region of 
maximum radial strain. A spherical bead thermistor 5 mils in diameter 
was attached to the underneath side of the dome electrode at the 
centerline with pressure sensitive tape. 

The strain gage and thermistor were connected into Wheatstone 
bridge circuits. Bridge excitation was one volt, in order to minimize 
self-heating of the devices. The thermistor was calibrated in an oil 
bath using a NBS-calibrated thermometer graduated in 0.01°C. The 
strain gage was calibrated in deflection of the dome by installing a 
small front-surface mirror (weight 0.204 gms) on a balsa wood plug 
in the dome. The transmitter was excited by the JRB circuit (a nearly 
constant current circuit) at various current levels and the steady- 
state deflection of the mirror was measured with an optical interfer- 
ometer. The strain gage bridge output was measured for each incre- 
ment of deflection. 

During each test, thermistor bridge output, strain gage bridge 
output, JRB circuit voltage across the transmitter, and current 
through the transmitter were recorded simultaneously on a four- 
channel Sanborn recorder. Results for four typical tests on a par- 
ticular transmitter are reproduced in Fig. 10. The test shown in lOd 
was a repeat of 10c changing onry the scale factor on the original 
recording equipment. 

In addition, detailed knowledge of the thermal resistance of granu- 
lar carbon in a brass container was necessary for the application of 
the analysis presented here. This thermal resistance was found to 
consist of two parts: an effective thermal conductivity for a con- 
tinuum representation of the granular carbon and a thermal contact 
resistance at the interface between the medium and the container 
walls. Experimental measurements of these two quantities have been 
recently completed and are reported in Ref. 4. 

In Fig. 10, two types of response are apparent. There is an early 
time transient for which the back electrode temperature can be as- 
sumed constant. After a few seconds the back electrode and the whole 
transmitter begins to heat up. The subsequent over-all growth coupled 
with the thermal expansion of the granular carbon itself combines 
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Figure 10. — Transient response of an excited Tl telephone transmitter, (a) Run 1, 
(I V) nv tra „siont = 0.056TF; (b) run 2, (IV)„ V transient = 0.156TF; (c) run 3, 
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in such a manner that the transmitter resistance (voltage) reaches 
a constant value. No attempt will be made here to describe analyti- 
cally this combination of effects at times later than about 6 seconds. 
However, the early time transient, which produces the resistance 
change, can be described. 

We first calculate N m for the standard Tl transmitter. Calculating 
the thermal conductance of the diaphragm from (36) 



= 5k,, cos i 12.5° = 66Q w/(m 2.o C) 
a In (so/s,) 



when 



5 = 0.0076 cm, 

k d = 140W/(m-°C), 

a = 0.28 cm, 

s =2.1 cm, 

8, = 0.286 cm. 

From Ref. 4 we take the thermal conductivity of the carbon to be 
0.24 W per (m-°C) and the thermal contact conductance as h c = 300 
W per (m 2 -°C). Consequently, the total conductance is the sum of the 
contact conductance in parallel with the conductance of the dia- 
phragm : 

lAtot = 1/660 + 1/300; 
hence, 

h lot = 206W/m 2 -°C). 
Using htot in N m we have 

N Bi = ha/k = 2.4. 

Turning to Fig. 6 we see that a${a, oo)//? = 0.19 and Fig. 9 gives t 
= 0.18 for this value of .V fli . Now, the thermal diffusivity k = k/ P c = 
0.32 X 10 -0 m 2 per second for granular carbon where P = 900 kg per 
m 3 and c = 840 J per (kg-°C) (for the loose state). Consequently, 

m r n (b - a) 2 = 14seconds> 

K 

where 6 = 0.44 cm from Fig. 5. In the determination of [3 we use the 
average values (IV) nv as obtained by a graphical integration of the 
transient portions of the appropriate curves of Fig. 10. Taking the 
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value from Fig. 10b as typical we have 

= (iv^ = 5l7 cm0C 

Therefore, = 0.19(3/a = 3.5°C. 

By comparing the time constant of 1.4 seconds to the measured 
values of t in the range of 3 seconds, one could say that this is rea- 
sonable agreement especially since the capacitance effects of the 
diaphragm and the electrodes have been neglected in the analysis. 
The capacitance in the real situation would tend to give higher values 
of t as observed. 

The calculated temperature rise of 3.5°C compared with the meas- 
ured "steady-state" value of 1.0°C needs some explanation. Of course, 
the capacitance and contact resistance of the thermistor would tend 
to give rise to the lower measured value. But more important is the 
fact that the predicted value is for the granular carbon continuum 
which is linked to the brass dome through a contact resistance. The 
temperature drop across that resistance can be calculated by multi- 
plying the temperature rise above ambient (6) by the ratio of the 
contact resistance (l/h c ) to the total resistance (l/h tot ) ; 



AT - aKl/n ' } = 2 4°C 

**■* continuum to thermistor ~ i n " -^ >'■ 



e(i/h e ) 



Therefore, LT mm ^, „bove ambient = 3.5°C - 2.4°C = 1.1°C; which is 
in remarkably good agreement with the measured value of 1.0°C. The 
other measured values, at the end of the early transient (namely, 
0.33°C, 1.56°C, and 1.61°C), scale linearly with the input power; 
hence, the agreement with analysis remains good for all the results 
given in Fig. 10. 

Consider briefly a transmitter made of invar. For invar the k d 
value is 18 times smaller than that of aluminum. Consequently, 

lAtot = 1/36.6 + 1/300; 

hence 

h M = 32.6W/(m 2 -°C) 

and N„i = 0.4. From Fig. 6 we see that the temperature rise would 
be 1.74 times that under corresponding conditions in the Tl trans- 
mitter. From Fig. 9 the value of t would be 30 percent longer. This is 
in qualitative agreement with the observations of Fischer and Gaudet. 2 
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Turning now to the displacement we first calculate the extension of 
the slant height of the aluminum diaphragm using (43). Taking v = 
1/3 and a t = 23.4 X lO-^O 1 we obtain for a temperature rise of 1°C 

u = -0.029 j«m. 

The minus sign indicates that the diaphragm grows toward the dome 
electrode. The expansion of the dome electrode is obtained from (37) 
using a t = 16.2 X lO-^C" 1 for brass. For d = 1°C, Ar = 0.04 /mi. The 
minimum displacement, given by (45) , becomes 0.0064 pm, while the 
maximum displacement given by (46) becomes 

M B „ = r 1 ^- - t^- = 0.180 + 0.134 = 0.32 /mi. 
tan a sm a 

The upper limit which assumes a perfectly rigid dome electrode com- 
pares reasonably well with the measured value of 0.18 /xm from Fig. 
10b. Notice that the top of the dome electrode moves an additional 
0.04 /mi. 

The most significant factor arising out of this rather crude model 
for displacements is that the thermal packing appears to be more 
sensitive to the thermal expansion coefficient for the dome electrode 
than any other part. Hence, a simple and inexpensive modification 
of the Tl transmitter would have been to change only the dome elec- 
trode to invar. Although the above analysis does not strictly apply 
to the modified transmitter design of Huffstutler, 1 the same conclu- 
sion with regards to an invar dome seems appropriate. 

The thermal expansion of the carbon produces an additional pack- 
ing effect. The change in volume caused by the heating could be 
calculated by using the expansion coefficient value 5 of 10.4 X 10" CO C _1 
and the local temperature as given in Fig. 7. An estimate of such a 
calculation yields a dome displacement equivalent to 0.05 fim in the 
direction opposite the displacements mentioned above for the condi- 
tions of Fig. 10b. 

V. CONCLUSIONS 

The resistance change of the Tl transmitter resulting from the 
joulean heating can be associated with the thermal expansion of the 
dome electrode and the aluminum diaphragm. The agreement of the 
model presented here with the measurements of Sandahl and Pretty- 
man attests to its validity for the early time transients. The displace- 
ments which produce thermal packing are a stronger function of the 
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thermal expansion coefficient of the dome electrode than any other 
part. A simple modification would be to make this part of invar in 
the Tl or any similar transmitter. The model also demonstrates that 
the displacements associated with the temperature rise of the dome 
electrode are affected less by changing the thermal parameters of the 
system than by reducing the thermal expansion coefficient in going 
from aluminum and brass to invar. The thermal contact resistance 
between the granular carbon and the dome electrode reduces the tem- 
perature rise significantly. 
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APPENDIX A 

Localized Heating in the Carbon Granules 

Because of the asperities of carbon granules and of the "bunching 
up" of the current flow at the contact areas, the importance of 
localized heating should be checked. It can be assumed that the 
localized effects are virtually all diffused at the values of time for 
which the Fourier modulus kU/D 2 is approximately unity. Then on 
setting the characteristic dimension of the granule, D equal to 0.2 
mm, and the thermal diffusivity, k, of solid carbon equal to 0.2 X 10" 4 
m 2 per second, the value of the time constant, t , becomes approxi- 
mately 

^o, ( 0.2) 2 X 10" 6 , . 
/o = Z)A= 0.28X10" 4 =1 - 4msec - 

We see that the thermal time constant for the solid carbon granules 
is at least three orders of magnitude less than the thermal time con- 
stant of the thermal effects being observed. Consequently, any lo- 
calized heating is rapidly diffused and the granular carbon can be 
treated as a continuum. 

APPENDIX B 

Transient Temperature Response of the Carbon Chamber 

Consider the transient portion of the heat conduction problem 
specified by (24) through (27) and repeated below for convenience. 
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d#/dr = d z &/dv 3 (50) 

*(,, 0) = F/0 ~ F ^ _ ( /_ F) J - v(v) (51) 

t?(0, t) = (52) 

d#/d v (l, r) 4- (1 "^ (N Bi + 1)0(1, r) = 0. (53) 

Assuming a separable solution of the form 

*(,,t) = X(r,)T(r) (54) 

so that 

X = A sin At7 + B cos Xtj (55) 

T = Ce" Xlr . (56) 

But (52) requires that B = and (53) specifies that 

\A cos X + ^A sin X = 
or, 

X cot X = -c, (57) 

where 

c, = (1 ~ F) (iV B ,- + 1). (58) 

To satisfy the nonhomogeneous initial condition we set 

*(l, 0) = Z 4„X n , (59) 

n = l 

where X„ = sin k»r) and X„ is the nth positive root of (57) . Multiply- 
ing both sides of (59) by X,„ , integrating, and using the orthogonality 
condition we have 

f XMdb, o) dt, 

A n = J -±—— — ■ (60) 

/ XI d v 

•'o 

Evaluating the integral in the denominator where the boundary con- 
ditions have been used we have 
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The transient part of the solution can thus be written as 

*0», t) = 2 r 2 2 ^ n t^ exp (-X'r) sin X„ij f tffo, 0) sin X n 77 dr,. 

(62) 

To evaluate the integral in (62) we first notice that (51) can be 
written as 

tJ(ij.O) = c, - c 2V + -z ° 3 

1 - C 4 77 

so that 

/ t?(rj, 0) sin X„77 ^77 = c /X n [l — cos X„] — c 2 /X„ ! [sin X„ — X n cos X„J 

"'0 

f Xn sin g (/£ 

+ c 'i„ x^a' (63) 

The last term which can be expressed in various forms is the one 
which causes the series to converge very slowly. 
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